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Abstract 

The variable phase approach to potential scattering with regular spheri- 
cally symmetric potentials satisfying (|l|), and studied by Calogero in his book 5 , 
is revisited, and we show directly that it gives the absolute definition of the 
phase-shifts, i.e. the one which defines 5n{k) as a continuous function of k 
for all k > 0, up to infinity, where <^(oo) = is automatically satisfied. This 
removes the usual ambiguity ±mr, n integer, attached to the definition of the 
phase-shifts through the partial wave scattering amplitudes obtained from the 
Lippmann-Schwinger integral equation, or via the phase of the Jost functions. 
It is then shown rigorously, and also on several examples, that this definition of 
the phase-shifts is very general, and applies as well to all potentials which have 
a strong repulsive singularity at the origin, for instance those which behave 
like gr~ m , g > 0, m > 2, etc. We also give an example of application to the 
low-energy behaviour of the S-wave scattering amplitude in two dimensions, 
which leads to an interesting result. 
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1 Introduction 



In quantum scattering theory with a spherically symmetric potential V(r), the 
phase-shift, defined for each partial-wave by the asymptotic behaviour of the radial 
wave function at large distances, has the usual ambiguity of ±7m, n integer. When 
the potential is regular, i.e. is L 1 (a, oo), and satisfies the Bargmann-Jost-Kohn 
condition 1 ' 2,3 

roo 

/ r\V(r)\dr < oo , (1) 
Jo 

one can show that the phase-shift Se(k) at infinite energy satisfies the condition 

tg Si(oo) = , or sin <^(oo) = , (2) 
so that one can make the "canonical" choice 



fc(oo) = , (3) 

and then proceed downwards by continuity for finite values of k. 

The same ambiguity exists, of course, when the phase-shift is defined through 
the phase of the Jost function 1-3 . The problem becomes even more serious when 
the potential is singular and repulsive near r = 0, and short-range otherwise : 



V(r) ~ g r m , g > 0, m > 2 ; ~ gexp(a/r), g > 0, a > ; etc. (4) 

r— > r^r 

In these cases, one can still define the S'-matrix, Si = exp(2i5e(k)), with a 
continuous phase-shift Se{k), which is again given by the asymptotic behaviour of 
the wave function for r — ► oo. However, now, one has 7-9 

5 e (oo) = -oo , (5) 

and this adds to the difficulty of finding a unique phase-shift. 

From the above remarks, it would seem therefore satisfactory to find a unique 
definition of the phase-shift itself by some formula containing V(r) and the wave 
function, and such that (|3]) is automatically satisfied when the potential satisfies (jl]). 
And then, to see whether such a definition can be extended, eventually with some 
modification, to potentials which are singular at r = or at r = oo (long range). 
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The answer to the above quest has already been found, although not written 
down explicitely in the form fl2D| ) we give below. It is given by the variable phase 
method 5 . To simplify the algebra, and see what is essential, let us take the case of 
the 5- wave, £ = 0. We have the radial Schrodinger equation 



<f"(k, r) + k 2 (p(k, r) = V(r) (p(h, r) 



r e [0,oo) , <p(0) = , y/(0) = 1 . 
If we write now the wave function as 



(6) 



(f = A(k, r) sin[/cr + 8(k, r)] 



A(k,r)>0 , A(k,0)^0 , 6(k,Q) = 
it can be shown that one has formally the differential equation 5 



' d 1 
— 5(k, r) = — — V(r) sm 2 [kr + S(k, r)] 
dr k 



5(k, 0) = 



and then 



A(k,r) 



k 



i J o r V(t) sin 2[kt+S(k,t)]dt 



(7) 



(8) 



(9) 



provided (§) has a solution, and the integral in (|9]) is finite. We shall see this in a 
moment. In essence, A is the amplitude of the wave function ip, whose oscillations 
are given by sin[/cr + 5(k,r)]. The function 5(k,r) can be interpreted as the local 
phase-shift since 8(k,R) is the phase-shift due to the cut potential V(r)8(R — r). 
The total phase-shift is defined as 



6(k) = lim 6(k,r) . (10) 

r— > oo 

In order to see whether (H) has a unique solution, we can write it as 

8(k, r) = — - [ r V(t) sin 2 [kt + 8(k, t)] dt . (11) 
k Jo 

One can then try to solve this nonlinear integral equation by iteration, starting from 
the zeroth order approximation 
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S^(k, 







(12) 



Numerical calculations show that this process is fast converging for usual regular 
potentials used in Nuclear Physics, and indeed has been used in practice. 
From (0), it is obvious that the phase-shift is given by 

I roc 

S(k) = -~ V(r) sin 2 [Ax + 6 (k,r)]dr . (13) 
k Jo 

This formula shows that, for short-range potentials, the tail of the integral : 



k Jr 



V(r) sin 2 [kr + S(k, r)] dr 



1 f°° 
< - / \V(r)\dr 
k Jr 



(14) 



can be made very small if we choose R larger than the "range" of the potential. 

In order to study the nonlinear integral equation ([11]) in a rigorous way, we can 
use the inequality 1,2 



sinxl < C 





X 




1 + 


x\ 



x real 



(15) 



where C is an appropriate constant. Noting that x/(l + x) is an increasing function 
of x for x positive, we have 



sin(x + y)\ < C 





x + y 




1 + 


x + 


y\ 



< c 





X 


+ 


y 




1 + 


X 


+ 


y\ 



(16) 



Using this in the integral equation (|TT|) , we find 



\5{k,r)\<^- f Q \V{t)\ 



kt+ \5(k,t)\ 



1 2 



dt 



k > . 



17) 



l + kt + \5(k,t)\ 

It is then obvious that an upper bound A(k,r) for \5(k, r)\ is obtained from the 
solution of the integral equation 



A(k,r) = ^£\V(t)\ 



kt + A(k,t) 



1 2 



dt 



l + kt + A(k,t) 

We can again solve this integral equation by iteration, starting from A^ ) = 0. It 
is obvious here that the solution cannot become infinite at any point on [0, oo). 
Indeed, if A(k,r) becomes infinite at r = r 1; then, as r f r 1; the fraction in ( |18| ) 
becomes one, and since V is supposed to be L 1 , we get a contradiction. 
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We show in Appendix A that the equation fll8|) has always a unique solution 
for all values of k > 0, provided V satisfies ([]]), and therefore that ( |TT| ) also has a 
unique solution. We show also that lim A(k, oo) = 0, from which we conclude then 

k— > oo 

directly that we have @. However, it is not really necessary to use the integral 
equation. Indeed, as we shall see in the next section, we can express the right-hand 
side of (|TT1) in terms of the regular solution ip given by (||), and its derivative <p', to 



obtain the new formulae 



w,r)=-*j;v® * , (i9) 



and 



*<*> = -*f vw **Tiv* ' (20) 

We shall see in the next section how to generalize these equations to higher t. 
Likewise, the amplitude A(k, r) can also be written as 5 

A(k,r) = - yV 2 + ^V . (21) 

rC 

Now, we know, quite generally, that for potentials satisfying ([]]), the wave function 
<p(k, r) exists for all values of k, real or complex, and all values of r > 1_4 . Likewise 
for (p'(k,r), and we have, of course 

<p'(k,0) = l , <p(k,r) = r + ro(l) . (22) 

r— » 

Moreover, we have also that ip and ip', for any fixed real k > 0, are real and con- 
tinuous function of r, and bounded as r — ► oo. Also, from the general theory of 
differential equations 12,13 , we know that, for any k, <p and <p' cannot vanish simulta- 
neously at any point r = > 0, since this would entail that (p = 0, a contradiction 
with (p'(k, 0) = 1. It follows that, because of ([!]), the integrals in (|I9|) and (p0|) are 
absolutely convergent for all real values of k > 0, and define continuous functions 
of k. Therefore, ( |T9| ) and define, in a very nice and simple way the "local" and 
the total phase-shift, respectively. We shall come in the third section to the case 
k = 0. 

Several remarks are now in order. First of all, (|20|) shows that, for a potential of 
a given sign, the phase-shift has the opposite sign, something well-known 1,2 , and also 
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obvious on ( |ITD and fll3l). Secondly, it is obvious on (|T9D and ( p0|) that the phase- 



shift is independent of the amplitude of tp. Multiplying tp by a constant factor X(k) 
independent of r leaves (|i9|) and ( p0|) invariant. This is as expected, of course. 
Finally, from fl2T|) we find 

A(A;,r)>0 , A(A;, 0) = t ^ (23) 

which was assumed in (|7]). Now, for r — > oo, the asymptotic behaviour of p and tp' 
are given by 1 " 4 

I FIT) I 

tp(k,r) ~ i-4^i sinffcr + 5) + o(l) , (24) 

oo 

and 



tp'(k,r) |F(A;)|cos(A;r + 5) + o(l) , (25) 

where -F(fc) is the Jost function, which never vanishes for k > 1_4 . Using now ( p4|) 
and (|D in (0), we find 

A(* )0 o) = !^jP , (26) 

and this, used in (0), leads again to (H), a consistency check. 

We have now to see whether, for regular potentials satisfying (ffl), the phase-shift 
defined by ( p0|) satisfies (||) : 5(oo) = 0. This is very easy to check. Indeed, for k 
real and going to infinity, we have uniformly in r 1,s ' 4 , for any r € [0, R], R < oo, 

sin kr 1 . „ , . 

^• r > ^ — + • (27) 

and 



<p'(fc,r) = cosfcr + o(l) . (28) 

oo 



Using these in (p0|), together with (0), leads to 



<5(A;) ~ —k [ V(r) S ^ n , ^ r dr = — [ V(r) - — C( ^ s ^ r 
fc^ oo Jo k 2 Jo 2k 



~ l Rv { r ) dr r sm2kt dt = - [ R sxa2kt dt ( [ 
Jo J o Jo \Jt 



R V(r)dr) . (29) 
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Now, in general, W(t) = j^Vij-) dr is L 1 (0, oo). Indeed 



\W(t)\dt= / dt 



oo 



V(r) dr 



oo 



< / dt \V(r)\dr 



roo rr poo 

= \V{r)\dr dt= r\V{r)\dr < oo (30) 
Jo Jo Jo 

by virtue of (11]). Therefore, in the last integral in ( p9|) we have a Fourier sine 
transform of an L 1 function. From the Riemann-Lebesgue lemma, it follows that, 
as k — > oo, it vanishes. Therefore, 5(k — > oo) = 0. This argument, which is quite 
general, can be made more precise (see Appendix B). It follows that the definition 
of the phase-shift given by ( p0|) in terms of the well-defined regular solution <p, (|6|), 
is indeed an absolute definition of 5(h) for all k > satisfying automatically (|3|). 

We shall see later that (|20|) can be extended to potentials which are outside 
the Bargmann-Jost-Kohn class, especially to those potentials which are strongly 
repulsive at r = 0. We are going now to give first the derivation of (f|), © and (^l|). 
These are found essentially in 5 , except for the integral forms ( |I9"D and (|20|), which 
are the basic equations of our paper. We reproduce the proofs for the convenience 
of the reader. 



2 Derivation of (19) and (20) 

We follow the usual method given in Calogero 5 for deriving the differential 
equation (|8|). Differentiating (^), we find 

iff = A' sin(kr + 5) + A(k + 5') cos{kr + 5) . (31) 

We have here two unknown functions A and 5, and only one equation, the Schrodinger 
equation (|6]), at our disposal. We can therefore impose a relation between A and 5. 
We impose 

A' sin(£x + S) + AS' cos(Ar + 5) = . (32) 
It follows that ( |3~TD becomes simply 

<p' = Ak cos(kr + 5) . (33) 
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Differentiating now this, and using the Schrodinger equation (^D, we find the new 
equation 

A'k cos(kr + 5) - Ak5' sin(£r + 5) = VA sin(kr + 5) . (34) 

Combining this with fl32|), we find the differential phase equation (|8|) and the am- 
plitude equation (^|), a complete set equivalent to the Schrodinger equation (]6|). 

However, we can go back again to ip and ip', using (0) and (|33|) . Adding their 
squares, we find 5 

A*(k,r) = l 2 (^ + kV) ■ (35) 
It follows that (H) can be written 



8' 



- l -V^\kr + 5) = -\ V^ 2 



-k V 



(36) 



ip' 2 + k 2 (p 2 

This is the basic equation, giving the "local" phase-shift 5(k, r) in terms of ip and 
ip'. Integrating it, we find 



6{k,r) 



<p 2 (k,t) 



~ k Jo V ^ cp' 2 (k,t) + k 2 <p 2 (k,t) 



dt 



(37) 



where the initial condition 5(k, 0) = has been used. Making now r — > oo, we 
get the total phase-shift, (p0|). We have already checked that in ([19]) and (|20| ) the 
integrals are absolutely convergent, and define continuous functions of k for all k > 0. 

Higher waves. 

The method for i ^ is quite similar 5 . We have to deal now with the radial 
Schrodinger equation 



ip"(k,r) + k 2 ipe{k,r) = 
and its free counterpart when V = : 



+ V{r) 



(38.a) 
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dr 2 \v t \vt. 



These free solutions are given by 



1-6 



ut(kr) 



ve(kr) 



(2£+l)\ 



lixkr 



'irkr 



J e+ i(kr) 



N £+ i(kr) 



1(1 + 1) fu? 
r 2 \v tj 



{kr) e+1 



T(2l + 2) 
2 e T{e+l) 



r- o (2£+ 1)!! 



(2£- 1)!! 



r-» (£x) £ 



+ ■ 



+ 



1 

> — 
2 



(38.b) 



(39.a) 



We have also 



Ui(kr) 
ve(kr) 

Their Wronkian is given by 



sin(/cr — ■^£' K ) ' 



cos(fcr £n) + ■ ■ 

2 



(39.b) 



W(ut, vi) = u' t V£ — ui v' e = k . (40) 
As for the regular solution ip e (k,r), it is customary to normalize it such that 1-4 



(pt{k,r) 



,£+1 



+ 



(41) 



r-o (2^ + 1)!! 

Now, again, under the condition (|T|) on the potential, one can show that both (p 
and ip' exist for all real or complex values of k, all r > 0, and are continuous in both 
variables 1-4 . Moreover, one has the asymptotic behaviours 



<p' £ (k,r) 



,=oo^P s^{kr- 1 -£n + S e (k))+o(l) 



\Fe(k)\ 
k e 



cos (kr - -£ir + 6t(k) \ + o(l) 



(42) 



where Fg(k) is the Jost function, well-defined and continuous for all k > 0. For k 
real, one has also, uniformly in r for any r G [0, R], R < oo, 



9 



<Pe(k,r) 



sin/cr o(l) 



oo k e+1 k i+1 
cos kr o(l) 



(43) 



k l 



We write now 



(fe(k,r) = A e (k,r)[u e (kr) cos 5 e (k,r) + v e (kr) sin 5 e (k,r)], (44) 

A £ (k,r) 7^ for all r > 0, <5^(A;, 0) = 0. We have now again two unknown functions 
Ag and 5i, and only one (differential) equation to determine them. We can therefore 
impose a relation between Ag and 5g. Simplifying the writing, we impose 



A'[u cos 5 + v sin 5] + A[— u sin 5 + f cos 5]<5' = 
Differentiating now (44) and taking into account (45), we find 

if' = A[u cos S + v ' sin S] 
One more differentiation gives us now 

ip" = A'[u' cos 5 + v' sin 5} + A[u" cos 5 + v" sin 5} 
+A[—u' sin 5 + v' cos 5] 5' . 
Using now (38. a) for <p", and (38. b) for u" and v", we find from (47) 

A'[u' cos 5 + v' sin 5] + A[—u' sin 5 + v' cos 5}5' 
= VA[u cos 5 + v sin 5] . 



(45) 



(46) 



(47) 



(48) 



We have now two equations, namely (45) and (48) to determine A and 5. We can 
write them, symbolically, as 



aA' + bA5' = 



(49) 



cA' + dA5' = aVA 
Eliminating A', and using A ^ and (40), we find 



(50) 



fO 



and 



ad — be = —(u'v — uv') = —k 



(51) 



S[ = — — V(r)[u£ cos 5 £ + i^sin^] 2 , 

K 

to which we have to add Si(k, 0) = 0. Likewise, eliminating 5', we find 



(52) 



kA' e = V(r)Ag[u£ cos 5g + vg sin Si] [—u e sin 8% + f^cos^ 



(53) 



As it is easily seen, we have again equations very similar to the equations for t = 0. 
Once (52) is solved, with the boundary condition Si{k, 0) = 0, we can replace its 
solution Si in (53), and integrate it to get Ai : 



A e (k,r) = A e (k, 0)exp 



1 

k Jo 



(54) 



However, as we saw previously, we can write both Se and Ai in terms of ip£. From 
(44) and (46), we can calculate sin<5£ and cos<5£ : 



sin S 



u'tf — wp' 



kA 



(55) 



cos S = 



Using now sin 2 S + cos 2 S = 1, we get 



vtp' — v'ip 
kA 



A 2 - — 

~ k 2 



(u e ip' e - u' e ip e f + (viip'i - v'iipiY 



Using also (55) and (56) in (52) and remembering (40), we get 



«=4"w-* r) 



k v 'Aj(k.r) 



Integrating this now, we have 



(56) 



(57) 



(58) 



6 e (k,r) = -k / V(t) 



(p 2 e (k,t) 



lo " vv [(u e (p' e - u' £ ipi) 2 + (vnp'i - V £ (pi) 2 } 
All these formulae, so far purely formal, are very similar to those for £ = 0. However, 
using now the behaviours of <p, ip', u, u', v, v' for r — > 0, we can check that all our 
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above formulae are meaningful, i.e. the integrals are convergent at r = 0. Again, it is 
clear from (57) that A^(k, r) 7^ for all r > 0. Indeed, if Ai is zero at r = r , we must 
have, according to (57), up' — u'(p = and vip' — v'(p = at this point. Calculating 
again <p and <p' from these two equations, we find that, if k 7^ 0, tp = <p' = at r = r$. 
And this entails that <p = everywhere, a contradiction with (41). Therefore, our 
assumption on A% is satisfied : At(k,r), for k > 0, never vanishes for r e [0, 00). 

Let us now look at the behaviour of 5e(k,r), (59), k > for r — ► 0, £ > —1/2. 
From (39. a) and (41), it is easily seen that the numerator behaves like r 2l+2 whereas 
the denominator becomes a positive constant. Therefore, we have, for k > 0, 

6 t (k,r) ~ q constant J\ 2i+2 V{t)dt = r 2l+l o(l) . (60) 

This shows that our assumption on 5e(k, r) is also satisfied : 8g(k, 0) = 0. Using (60) 
in (53), we find also that we have a convergent integral in (54). This completes the 
validity of the method and its consistency. 

We have now to look at (57) and (59) for r — ► 00. From (39. b) and (42), we find 
easily, for all k > 0, that we have 

where Fe(k) is the Jost function, finite and continuous for all k > 1_4 . From (59) 
for r — > 00, we find 

fc(jfe) = -Jfc rV(r)— n!^'^, tt^t rfr . (62) 

Jo v ; [(u' t <pt - u e ip' e ) 2 + (v' e pe - v e <p'e) 2 ] 

and it is easily checked from (39. b) and (42) that the integral here is well-defined and 

absolutely convergent, and since (p, ip', u, u', v and v' are all continuous functions 

of k for all k > 1_6 , the same is true for 5t(k) : the phase-shift is a continuous 

function of k for all k > 0. To check (3), we can use (43) in (62), and we find, as 

for the case i = 0, that 5e(oo) = 0. Moreover, k e+1 Ae and 5e are also continuous 

functions of k for all k > 0. They are given, respectively, by (57) and (59). The 

physical phase-shift is given by (62), an absolutely convergent integral for all k > 

under the assumption (1), and one has also (3). Therefore, (62) gives an absolute 

definition of the phase-shift 5i(k) in general for all k, and for potentials satisfying 

d]]). We shall see in the next section that (|20|) and (62) are valid also for singular 
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repulsive potentials as well. The case k — will be considered at the end of the 
next section. 

Remark. Making I = 0, one finds, as expected, all the formulae we found previously 
for the S'-wave. Before ending this section, let us mention that by combining (38. a) 
and (38. b), and integrating from to r, using the appropriate boundary conditions 
at r = 0, namely (39. a), (40), and (41), to evaluate the integrated terms, one can 
calculate u'tp — mp' and v'ip — v<p f , so that (62) can be written, for £ > — \, also as 

W) = -k T Vir)-^ Z^TZ r T^dr . (63) 

Jo (J u t (ptVdt) 2 + {k f + / !)< (feVdt) 2 

Note here that we have now well precised boundary conditions at r = for all 
the functions entering in (63), so that care must be taken in using it, whereas in 
(20) and (62), the normalization of (p at r — does not matter. For — | < t < 0, 
as we shall see in the next section, (62) and (63) are valid provided we use there 
the "distinguished" pure Bessel solution for <p, given by the integral equation (68). 
Indeed, now, both free solutions ii£ = \^rJ e+ i(kr) and ve = \/rN t+ i{kr) vanish at 
r = 0, and so, as free solution, we can start from any combination au£ + bvg, and 
use it as the inhomogeneous term in (68). We get then always a solution ipg with 
(p(k, 0) = 0. The "distinguished" solution is the one with 6 = 0. 

3 Domain of validity of (20) and (62) 



As we have seen, formulae (20) for the S'-wave, or its generalization (62) for 
higher waves, are valid for all k > and all i > 0, provided the potential satisfies 
the integrability condition (1) : rV(r) G L 1 (0, oo). Roughly speaking, this means 
that V is less singular than r~ 2 at the origin. We may now ask whether they are 
also valid for potentials having stronger singularity there, for instance V(r) ~ gr~ m , 
m > 2, g > 0, or gexp(^), g > 0, a > 0, n > 0, etc., as r — > 0. Rather than 
developing the general formalism for such general singular potentials (singular and 
repulsive at the origin), we shall consider explicit examples, and leave the full theory 
for a forthcoming paper. 
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i) We consider the formula (20) for the S'-wave, and take boldly the potential to 
be the centrifugal barrier 

V(r) = e -^±, £>0, (64) 

in the Schrodinger equation (6). The wave function is now just, up to an unimportant 
constant multicative factor, p = \^krj£ + i/ 2 (kr) which is of the form 4>(kr). It is 
easily seen that (20) is well-defined because the integral is convergent (absolutely) 
both at r = and r = oo. If we make the change of variable z = kr in it, we find, 
by writing <// = k-^cj) = k<j>(z), ' = d/dz, 

5= -£(£+!) [°° -. — ^ dz , (65) 

where = y/zJi + i/2(z). This last formula is now independent of k. The same 
is therefore true for the original formula (20) with our p(k,r). We can therefore 
calculate it at any value of k, for instance for k — 0. Using (39. a), we find 

s =^ t+ »rw*wT*\ = - t r£e = - t '* ■ (66) 

This is exactly the phase-shift of the centrifugal barrier £(£ + l)/r 2 since, without 
this potential, the wave-function is sin kr, and with the potential, <p = constant x 
sin(/cr — \£ir), as r — > oo, according to (42). 

In conclusion, our formula (20) is valid for repulsive singular potentials g/r 2 , 
g > 0, which violate (1) both at r = and r = oo. 

ii) We consider now the previous example, but with —1/2 < £ < 0. Here, 
proceeding as before, we find, as expected, again (66). Note that 5 is now positi- 
ve because the potential £(£ + l)/r 2 is negative. For £ = —\, the full solution is 
ip = \/krJ (kr), so that ip ~ ^fr and ip' ~ l/(2y / r) as r — > 0. 

Remark. Formula (20) was proved for regular potentials. In all rigor, in order to 
apply it to the centrifugal barrier potential £(£+ l)/r 2 , we must first regularize this 
at the origin, for instance by cutting it by 0(r — e), and then make e [ 0, or by 
replacing r in the denominator by (r + e), and again take the limit £ J, 0. However, 
as we saw, at the limit, we have already an absolutely convergent integral for all 
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£ > — |. Also, the derivation of (20) was based on the assumption (p'(0) = 1, i.e. 

(p(r) ~ r + • • • for r — * 0. However, for £ 7^ 0, we have rather ~ r e+1 -\ . The 

extra factor r £ comes from the regularized formula for the amplitude A in the limit 
e I 0, as can easily be seen 5 . 

iii) We consider now (13) with 

V(r) = y£±l± + V 1 (r), t>-\ , (67) 

assuming that rV\{r) £ L 1 (0, 00). Consider first I > — \. As we saw, the presence of 
this "weak perturbation", as compared to £(£+!) /r 2 , does not modify the behaviour 
of the regular solution ipe at r — 0, given by (41), to be compared with (39. a) 1-4 . 
And since (20), as we just saw, works with £(£ + l)/r 2 , £ > —1/2, applied to the 
regular solution <pe, it should work also when we add Vi, provided we use always pure 
Bessel functions as free solutions. This means that the Volterra integral equation 
which combines the Schrodinger equation and the boundary condition at r = is 

<Pe(k, r) = ^u e (kr) + J r G e {k; r, r') V(r') <p t (k, r') dr' , 

(68) 

G e (k;r,r') = \ [u t (kr)vt(kr') - u e (kr')v e (kr)] . 

We consider now I — — \. As has been shown in reference 10 , in order to formulate 
a decent scattering theory leading to the asymptotic form (42) with a well-defined 
Jost function Fe(k) and a well-defined phase-shift 5^ one has to make stronger 
assumptions on V than (1), namely 

/ r\V(r)\(l + \\ogr\)dr < 00 , 
Jo 

(69) 

PCX) 

/ r\V(r)\(\ogr) 2 dr < 00 , a>0 . 

i J a 

This will be used later for £ = —1/2, and k — > 0, where we give more details. 

In conclusion, formula (20) is valid for (67) and the regular solution There- 
fore, we have now two methods to deal with (67). The first one is to apply (62) to 
Vi, and the second one to apply (20) to the full potential (67). In any case, we get, 
of course, the full phase-shift 



15 



5f tal = 5,- l -h , (70) 

5e being the physical phase-shift due to V\, which is what interests us in scattering 
theory. Remember that, in (20) or (62), <p> is always the full solution : solution of 
(38. a) with V = V±, or (6) with (67), always together with (41). 

iv) We consider now more singular potentials, namely V(r) = g/r m , m > 2, 
g > 0. Here, it is known that the solution of the Schrodinger equation (1) which 
vanishes at r = behaves there as 7-9 

ip e (k,r) ^ o [\/(r)]- 1 / 4 exp(- f[V(t)]^d?j 
= Mr)=Q- 1/4 r m / 4 exp 



i— 2 I 

■Vfi 



2 m — 2 



(71) 



independent of k and i. The wave function (p and all of its derivatives y?^ vanish 
exponentially at r = 0. Notice that we can omit the factor g~ l l 4 in front of the last 
expression since our formulae for the phase-shifts are homogeneous in (p. In fact, 
at k = 0, the Schrodinger equation is soluble exactly, and its solution is 7 , up to an 
unimportant constant multiplicative factor, 

<p e (h = 0,r) = V^K^ (ML r -^\ , (72) 
m-2 — 2 J 

where K v is the modified Hankel function 6 . Using the asymptotic behaviour of 
K v (x) for x — > oo, we find indeed, up to constant multiplicative factors, the be- 
haviour shown in (71). On the other hand, we know that now, because of the strong 
singularity of V at r = 0, the phase-shifts does not go to zero as k — * oo, contrary 
to the case of regular potentials satisfying (1). Rather, one has the high energy 
behaviour 8 ' 9 



where 



6 e {k) = -A g™ k 12 ^ 1 H (73) 

oo 



2 r(3/2-l/m) • 1 ] 
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Since the main term in (73) is independent of £, we shall consider the case £ = 
to simplify the algebra, and therefore use (20). In this formula, the integral can be 
split into /(f + Jr. For the second integral, we have 

i*r-Hf£-°(^*) • ™ 

e very small (> 0), as k — > oo, provided i? = 0(A;~™ +£ ). The contribution of (75) 
can therefore be neglected if we compare with (73). Note that R — > as k — > oo. 
In the first integral, we can therefore replace, in first approximation, </? by O given 
in (71), and independent of k. We find then 




Making now the change of variable x = k 2 ^ m r, letting k — > oo, and noting that 
k 2 / m R — > oo, we find 

/ • • • ~ -<?£; m / — — - = —g ' k ™ / . 77 

Jo oo Jo g + x m Jo 1 +t m 

It follows that the approximate value of behaves asymptotically as 

5 app (k) ~ -Bg^k*^ , S = r . (78) 

fc^ oo Jo 1 + t 

This coincides with (73) up to a numerical factor, and is obtained without much 
effort, as we see. 

The above argument to obtain (78) is, of course, heuristic because (71) is uni- 
form in k only for k G [0, K], K < oo. But it can be made more precise and quite 
rigorous. In a forth-coming paper, we shall develop a technique to deal with all 
these problems in a unified way. 

Remark. The remark at the end of ii) applies here too. We must regularize first 
the potential, and then let £ J, 0. The exponential decrease of ip as r — > comes 
then from the amplitude A in the limit e j 0. 
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A low-energy example. 



Our last example is the low energy behaviour of the phase-shift for i = —1/2, 
where the potential is assumed to be repulsive (V > 0), and to satisfy the integra- 
bility conditions 

/ r|V(r)|(l + | logrlWr < oo , (79) 
Jo 



and 



roo 

/ r\ V(r) | (log rfdr < oo , a>0 . (80) 

J a 

This corresponds to the S'-wave Schrodinger equation in two space dimensions, and 
is interesting to study 10 . 

As we saw in iii), we must use here only the pure "Bessel" solutions. This means 
that the solution (p is the solution of the integral equation (68) for £ = — \ : 



(p(k,r) = Vkr J (kr) + / G(k,r, r') V(r') (p(k,r r ) dr' , (81) 

Jo 

where the Green's function is given by 

G = ^V^? [Jo(kr) N (kr') - J (kr') N (kr)\ . (82) 

The occurrence of logr in (79) is due to the presence of log (At) in N . Here, ip is 
normalized somewhat differently, but we know that the normalization of f does not 
matter in (20), (62). In (63), (f is normalized now as to ip ~ y/r for r — > 0. 

It is then shown in the above reference that the phase-shift, i.e. the phase- 
shift So(k) of the S'-wave in the two dimentional space problem, has the universal 
behaviour 

5 (k) ~ — + , (83) 
K ' o | log A; | 

i.e. the main term is independent of the potential. 

We are going to find (83) by using the formula (63) for i = —1/2, in which (p e is 
given as above, and, for z — > 0, 6 



18 



Uz) = l-- + 



(84) 



N (z) = - 

7T 



log - + 7 



Jn(z) +0(Z 1 



From the integral equation for ip, it can be easily shown that, under the assumptions 
(79) and (80), the low-energy behaviour of ip is given by <p ~ VkrJ (kr) ~ V^r + 

k— + 

• • •, which we have to normalize to ip ~ y 7 ?, and, for w and t>, in order to comply 
with (40), we must take 



•u ~ \fkr- + • • • 
t> ~ — v^At log fcr + • • • 

Using now the above low-energy behaviours in (63), we find 



(85) 



S(k) ^ 





rV(r) 



o (/ r tV(t)dtf + (1 - Jo tV(t) log kt dtf 



dr . 



(86) 



Since V was assumed to be positive, we can introduce the new variable X = X(r) 
by 



X(r) = \ log k\ T tV(t)dt , dX = \ \ogk\rVdr , 
Jo 



(87) 



which is a one-to-one mapping from r G [0, oo) to X G [0, |log/c|A], where A 
r\V(r)\dr. Letting now k — > 0, we find easily from (86) 



fc^ o 



1 



log A; | 



(88) 



Here again, as for the case of singular potentials, we do not get the exact constant 
Ti / 2 because we use ip ~ y/r, which is not uniform on the entire r-axis. We shall 
come back to this in more detail in the forthcoming paper. 

The case k = 0. 



Let us consider now k = in (20). It is known that, under (1), one has 



1-4 



<p(0,r) = y?'(0,oo)r + D + o(l) , 



(89) 
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where <p'(0, oo) is finite or zero, and D is finite also. Making now k = in (20), we 
find, for the S-wave scattering length 1 ' 2 

fc- o fc io y? 2 (0,r) V ; V ; 

At the origin r = 0, there is no convergence problem since <p'(k,0) = 1 for all k. 
However, at r = oo, because of (89), we must assume 

roo 

/ r 2 \V{r)\dr < oo (91) 
Jr 

in order to secure proper convergence, which is also well-known 1 ' 2 . But this is not 
yet the end. We must also be sure that <//(0, r) does not vanish for r > 0. This is 
surely the case if V is positive 1 ' 2 , but cannot be guaranteed otherwise. In conclusion, 
(91) is valid only when <//(0, r) ^ for all r > 0. For higher waves, in order to have 
proper convergence at r = oo, one needs 1 ' 2 

/ r 2e+2 \V{r)\dr < oo (92) 
Jo 

and again the non- vanishing of the denominators in (20) or (62). We shall see in a 
forthcoming paper, how to modify these formulae in the presence of bound states. 
However, we know that in all cases, Se(k) is continuous down to k + 0, and one has 
the Levinson theorem 5(+0) = nir, where n is the number of bound states 1-5 . 

Two-potential case. 

Formulae (62) and (63), as it is obvious, can of course be applied in the case 
where we have two potentials : 

V = V 1 + V 2 , (93) 

both satisfying (1). Here, we can apply either (20) to V, or (62) and (63) to 
where ut and vt are now replaced by two appropriate independent solutions of the 
Schrodinger equation ipi and tpi, normalized according to (40) : 

W[i Pl ,^ 1 ] = ^ 1 -^[ = k . (94) 
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With (20), we would get the full phase-shift, and with (62) or (63), the phase-shift 
8 2 due to V 2 . 
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Appendix A 



We have to study here the integral equation (18) 



A(*,r) = £j£V<t)| 
We solve it by iteration, starting from 



kt + A(k,t) 
l + kt + A(k,t) 



dt 



(A.l) 



A^\k,r) = . 



(A.2) 



Since x/(l + x) is an increasing function of x for x > 0, we get the increasing 
sequence of iterations 



where 



A w (k,r) < A (2 \k,r) < ••• , 



(A.3) 



A(1) ^ = ?/ >Ki^) * • 



AW(M = f/V(0l(y 



kt + A^-^M) ' 

+ + A( n " 1 )(A;,t) / 



(A.4) 



i) Assume now first that V is integrable at r = and therefore is L 1 (0, oo). It 
is then obvious that the increasing sequence of iterations is bounded by 



A(k,r) = ^ r \V(t)\dt , r>0,A;> 
k Jo 

It has therefore a limit, and one has the solution 







(A.5) 



A(k,r)= lim A {n \k,r) < A(k,r) 



(A.6a) 



for all r > 0, and all k > 0. For r — > oo, the same statement is valid. Indeed, if we 
note that both A(k,r) and A(k,r) are increasing functions of r, it follows that we 
have, when r — > oo, 



A(fc) = r lim o A(A;,r) < A(fc) = r fim o A(A;,r) = j J™ \V(t)\dt . (A.6b) 
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It is then obvious from (A. 5) that we have A(k — > oo) = 0, which gives in turn (3), 
as expected. 

ii) If we have only (1) : rV(r) G L 1 , we must refine slightly our argument. Since 
our problem is now the convergence of the integral in A^ n \k, r) at r = 0, we consider 
r very small. Now, as it is obvious, another sequence of upper bounds for A( n )(/c,r) 
is obtained by using the drastic inequality x/(l+x) < x in (A.l). We obtain in this 
way a sequence of upper bounds for A( n ) from the sequence of iterations of 

A(k,r) = ^ [ r \V(t)\ \kt + A(k,t)] 2 dt . (A.7) 
k Jo L J 

Putting A = kcu, we get 



u(k,r) = D \V(t)\(t + LuYdt , (A.8) 
Jo 

which is in fact independent of k, and must be iterated now. The above equation is 
nothing else than the integral equation for minus the local scattering length 

a(r)= lim^i (A.9) 

for the potential —D\V(r)\, and has been studied in the book of Calogero 5 , chapters 
11 and 12, where it is shown that the iteration of (A.8) leads to an absolutely 
convergent series expansion for the solution, provided r is small enough and rV(r) G 
L 1 . Alternatively, (A.8) is nothing else than the Riccati equation 

u'(r) = D\V(r)\[r + u(r)] 2 (A.10) 

with lj(0) = 0, which has been also thoroughly studied in the books of Hille 12 and 
Coddington and Levinson 13 , to which we refer the reader, with the same conclusion. 

Once we secure the solution of (A.7) in a small interval [0,r ], with A(k,0) = 0, 
we can then start at r = r , consider instead of (A.l) the integral equation 



A(k,r) = A(k,r ) + jJjV(t)\ 



kt + A(k,t) 



2 



dt , (A.ll) 
l + kt + A(k,t)\ y ' 

and proceed as before, by iteration. Here, we need only V G -^ 1 (ro, oo). We get then 
again an increasing sequence bounded by 

A(k,r ) + jf\V(t)\dt , (A.12) 
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with the same conclusions as before for the existence of the limit of the iterations, 
... etc. However, the high energy behaviour of S(k,r) cannot be obtained from the 
above analysis. 

A different method, which superseeds the above considerations, and provides at 
the same time the high energy limit of S(k, r) and of 5(k) = 5(k, oo) is as follows. It 
consists in neglecting A in the denominator of the right-hand side of (A.l) : 



A(k,r)<j£\V(t)\ 



kt + A(k,t) 
1 + kt 



D r r r 



kt + A 
1 + kt 



dt . (A. 13) 



Writing now A = ru, the two ends of (A. 13) lead us to the integral equation 



D 



CO 



(k + cuf 



dt 



k Jo ^"(i + kty ~" (A ' 14) 

whose solution provides still a stronger upper bound for A{k, r). Now, it is trivial to 
solve (A. 14). We just differentiate it, and integrate the differential equation, taking 
into account u(k, 0) = 0. This solution is just : 



u(k, r) = 



k I(k,r) 
1-I(k,r) ' 

t\V(t)\ 



(A.15) 



/o (i + kty 

The solutions exists as long as I(k, r) < 1, that is, in the interval [0, r ), where r (k) 
is given by 



D 



pro r 

Jo (1 



r\V(r) 



dr = 1 



(1 + kr) 2 

Note that, in any fixed interval [0,R], R < oo, we have 



(A.16) 



lim I(k,r) = 



oo 



(A.17) 



Indeed, from (A.15), and using 1 + kt > 1 and 1 + kt > kt, we have 



I(k,r) < I(k,oo) < D 



t\v(t)\ 

o 1 + kt 



dt 



re j roo 

= J(k, oo) < / t\V(t)\dt+-D \V(t)\dt , 

Jo k Je 

(A.18) 
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where, in obvious notations, J(k,r) is defined by D J r [t|V A (t)|/(l + kt)\dt. Now, we 
can make first e small enough, independent of k, in order to make the first integral 
in the right-hand side as small as we wish. Once e is fixed, we can than make k 
large enough in order to make also the second integral as small as we wish. This 
proves (A. 17). 

Remark. It is obvious that (A. 17) is uniform in any finite interval < r < R. 



As a consequence of (A. 17), r (k) defined by (A. 16) satisfies 



lim r (k) = oo , 



k—t oo 



(A.19) 



so that, the larger k is, the larger is the domain of validity of (A. 15). In any case, by 
making k large enough, and combining (A. 15) and (A. 17), we can have, uniformly 
in r in [0, R], 



u(k,r) <2k I(k,r) , k>K , 



(A.20) 



and we know that ru(k,r) is an upper bound for A(k,r). We can therefore, for k 
large enough, use (A.20) in (A.l) : 



kt + 2kt I(k,t) 



dt 



1 + H - = J(k, r) + 2D I(k, t)dt < J(k, r) + 2I(k, r) J(k, r) . 

(A.21) 

Now, as we saw for (A. 17), both / and J go to zero as k — > oo, uniformly in r for r 
in any finite interval [0, R] . Therefore 



lim A(k,r) = , < r < R . 



k^r OO 



(A.22) 



Once we have shown the existence of the solution A(k, r) in [0, R\, we can proceed 
as for (A. 11), write, for r > R, 



A(k,r) = A(k,R) + jf R \V(t)\ 



kt + A(k,t) 
l + kt + A(k,t) 



dt 



(A.23) 
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and proceed again by iteration. This way of proceeding is legitimate since (A.l), or 
(18), was obtained from the differential equation (8) and the bound (16). We can 
therefore start at any point r = R, and integrate after, provided we know 5(k,R). 
And we obtain a bigger upper bound if we replace \5(k, R) \ by A(k, R). The process 
is now very similar to what we had in i) . We get an increasing sequence of iterations 
A( n )(£;,r), with A^(k,r) = A(k,R), and a global upper bound for all A^ n \k,r) : 

A (n \k,r)<A (k) = A(k,R) + ^r\V(t)\dt , (A.24) 

k jr 

valid for all r > R, including r = oo, and all n. The sequence has therefore a limit, 
and this limit provides the solution A(k,r) of (A.l) for all r > R, including r = oo, 
with R < r (k), and r (k) given by (A. 16). Since we had also proved the existence 
of the solution in [0, ro], we have therefore proved the existence of the solution of 
(A.l) for all r > 0. Obviously, we have 

A(fc,0) = . (A.25) 

Therefore, we have, for r > R, 

A(k,r)<A(k,R) + ^r\V(t)\dt . (A.25) 

k Jr 

If we note also that A(k, r) is an increasing function of r, we secure the existence of 
A(k) = A(k, oo), and from (A. 22) and (A.24), we obviously have, for the phase-shift 



lim 5(h) < lim A(k) = lim A(k, oo) = . (A.26) 

oo oo k— > oo 

This completes our proof for the general case where rV(r) e -^ 1 (0, oo). 
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Appendix B 



Here, we shall show that, for regular potentials satisfying (1), one can make (3) 
more precise if one knows the detail behaviour of the potential when r — > 0. One 
has, indeed, the following : 

Theorem 1. Let us assume that V is continuous and bounded, away from the 
origin, in [0, R], R < oo, and is such that 



lim r 1+a V(r) = V , < a < 1 . (B.l) 



We have then, as k — > oo, 

m = Z^ cm (^y {1 _ a) _}_ + ... , (B.2) 

and conversely : (B.2) entails (B.l). 

We shall give the proof for the S'-wave, I = 0. The proof for higher waves is 
quite similar. Because of (14), we can just limit ourselves to (29) : 

5(h) = - J*dt sin 2kt (j* V(r) dr^j . (B.3) 

For the asymptotic behaviour of this integral (as k — > oo), we can now use the 
following theorem of Titchmarsh 11 : 



Theorem 2. Let f(x) and f'(x) be integrable over any finite interval not ending 
at x — ; let x 1+a f'(x) be bounded for all x, and let f(x) ~ x~ a as x — > 0. Then, 
denoting by F the limit of x 1+a f'(x) as x — > 0, we have 



00 7TQ; 1 

/(x)sinA;a; = F T(l - a) cos — — - (1 + o(l)) , (B.4) 

as fc — > 00. The converse theorem is also true if we deal with finite intervals in (B.4) 
since we can define the Fourier inverse transforms in a straightforward manner. This 
is indeed the case in (B.3). 

Our theorem 1 follows now immediately from the theorem of Titchmarsh applied 
to (B.3). Moreover, it is a Tauberian kind theorem, i.e. its converse is also true 
if we remember that the Fourier transform of an L 1 function f(x) is a continuous 
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function F(k) of k. More refined theorems containing logarithmic terms can also be 
proved. Examples are treated in Calogero's book 5 . For higher waves, the proof is 
similar by using (39. b) and (43) in (62). 



28 



References 



1. R. G. Newton, Scattering Theory of Waves and Particles, 2nd ed. (Springer- 
Verlag, Berlin, 1982). 

2. A. Galindo and P. Pascual, Quantum Mechanics, 2 volumes ( Springer- Verlag, 
Berlin, 1990). 

3. K. Chadan and P. Sabatier, Inverse Problems in Quantum Scattering Theory, 
2nd ed. (Springer- Verlag, Berlin, New York 1989). 

4. K. Chadan in K. Chadan, D. Colton, L. Paivarinta and W. Rundell, An Intro- 
duction to Inverse Scattering and Inverse Spectral Problems (SIAM, Philadel- 
phia 1997). 

5. F. Calogero, Variable Phase Approach to Potential Scattering (Academic Press, 
New York 1967). See especially chapter 3, and chapter 6, formulae (3) and 
(8). A complete list of references is given in this book. 

6. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions. 

7. W. Frank and D. J. Land, Singular Potentials, Rev. Mod. Phys. 43, 36-98 
(1971). 

8. N. Froman and K. E. Thylwe, Setting the question of the high energy behaviour 
of phase-shifts • • • , J. Math. Phys. 20, 1716-1719 (1979). 

9. O. Brander, High-energy behaviour of phase-shifts for scattering from singular 
potentials, J. Math. Phys. 22, 1229-1235 (1981). 

10. K. Chadan, N. N. Khuri, A. Martin and T. T. Wu, Universality of low-energy 
scattering in 2 + 1 dimensions, Phys. Rev. D58, 025014-1-19 (1998). 

11. E. C. Titchmarsh, Theory of Fourier integrals, 2nd ed. (Oxford University 
Press, Oxford 1962), pp 172-173, theorems 126 and 127 combined together. 

12. E. Hille, Lectures on Ordinary Differential Equations (Addison- Wesley, Read- 
ing, MA, 1969). 



29 



13. E.-A. Coddington and N. Levinson, Theory of Ordinary Differential Equations 
(McGraw-Hill, New York, 1955). 



30 



